Low temperature transition to a superconducting phase in boron-doped silicon films 

grown on (OOl)-oriented silicon wafers 
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We report on a detailed analysis of the superconducting properties of boron-doped silicon films 
grown along the 001 direction by Gas Immersion baser Doping. The doping concentration cs has 
been varied up to ~ 10 at.% by increasing the number of laser shots to 500. No superconductivity 
could be observed down to 40mK for doping level below ~ 2.5 at.%. The critical temperature T c 
then increased steeply to reach ~ 0.6K for eg ~ 8 at%. No hysteresis was found for the transitions 
in magnetic field, which is characteristic of a type II superconductor. The corresponding upper 
critical field H C 2(0) was on the order of 1000 G, much smaller than the value previously reported 
by Bustarret et al. [T] 



PACS numbers: 73.61. Cw, 72.62.Dh 

The discovery of a superconducting transition around 
40 K in the MgB2 compound [2] revived the interest for a 
specific class of superconducting materials : the covalent 
metals [3J 0]. Indeed, the high T c value is here a direct 
consequence of the strong coupling of the boron cr— bands 
(shifted up to the Fermi level due to the presence of Mg 2+ 
ions) with phonons. This system was the precursor of 
new covalent superconductors among which boron doped 
diamond C:B |5J, silicon Si:B pQ, silicon carbide SiC:B [BJ, 
as well as gallium doped germanium [7] . 

Eventhough these systems share the similarity of being 
heavily doped semiconductors for which superconductiv- 
ity appears in the partially unfilled valence band, they 
also present striking differences. In diamond, the onset 
of doping-induced superconductivity coincides with the 
metal-insulator transition (MIT) [3]. On the other hand, 
although silicon becomes metallic for boron doping on 
the order of 80 p. p.m., a concentration of several per 
cent was necessary to trigger superconductivity pQ . Such 
a concentration exceeds the solubility limit and can only 
be reached using sophisticated out-of-equilibrium doping 
techniques. This explains why a systematic study of the 
evolution of T c vs boron doping in Silicon is still lacking. 

We present here such a study for a series of high quality 
boron-doped silicon films. No superconductivity could 
be observed down to 40mK for doping levels up to ~ 2.5 
at.%. T c then rises sharply with boron concentration 
reaching ~ 0.6K for Cb ~ 8 at.%, a twice larger tem- 
perature than reported previously pQ. A similar rapid 
increase in T c vs doping has been reported in B doped di- 
amond films for which T c scales as (cs/c c -l) [8], where 
c c is the critical concentration corresponding to both the 



MIT and the onset of superconductivity. As shown be- 
low a similar dependence may also roughly apply to the 
present results but would require the use of a c c value 
on the order of 2.5 at.%, i.e. well above the MIT dop- 
ing threshold. We also show that, despite their higher T c 
values, the upper critical field H C 2(0) of our high quality 
films is about four times lower than the one previously 
reported in pQ. 

In order to overcome the solubility limit, a set of sam- 
ples have been prepared by Gas Immersion Laser Doping 
(GILD) [TJ A precursor gas (BCI3) is chemisorbed 
on a (OOl)-oriented silicon wafer which is subsequently 
melted using an ultraviolet laser pulse. During each 
melting/solidification cycle, boron diffuses from the sur- 
face into molten silicon, and is incorporated at substi- 
tutional sites of the crystal as the liquid/solid interface 
moves back to the surface of the epilayer upon cooling 
after the laser pulse. The in-plane lattice match to the 
substrate and the local incorporation of boron atoms of 
lower covalent radius result in the formation of biaxially 
strained pseudo-morphic epilayers. The B concentration 
was progressively increased by increasing the number of 
laser shots, from 1.4 at.% (50 shots) up to ~10 at.% (500 
shots), the laser energy being regulated in order to keep 
the doped layer thickness constant during the whole pro- 
cess, here on the order of 80-90 nm. Since our previous 
publication pQ, new laser optics have strongly improved 
the spatial homogeneity in energy within the spot and 
an ultrahigh residual vacuum (p ~ 10 -7 Pa) is reached 
in the reaction chamber. 

Measurements of boron concentration depth profiles 
were done by Secondary Ion Mass Spectrometry (SIMS) 
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FIG. 1: Boron atomic concentration depth profiles deduced 
from SIMS measurements for the indicated numbers of lasers 
shots in the GILD process. 



using a CAMECA IMS4f instrument. Reliable quantifi- 
cation of boron in silicon for such high concentration lev- 
els above 1 at.% is a major challenge for SIMS analysis 
due to matrix effects |10j as complex ionization and/or 
sputtering mechanisms may occur during secondary ion 
emission. We have shown that these effects remain small 
with an O^" beam in presence of a saturated oxygen 
flow at the surface [TT]. A primary 4keV O^sputtering 
beam of 90nA is rastered at oblique incidence (46.8°) on 
250x250 /im 2 of the surface. To avoid crater edge ef- 
fects, 10 B + , n B + , and 29 Si + ions coming from a central 
part of 10/xm in diameter have been selected and ana- 
lyzed simultaneously by an original Isotopic Comparative 
Method (ICM) (for details see [12]). Atomic concentra- 
tions are presented in Fig. 1 as a function of the depth 
measured by mechanical profilometry with a Tencor P10. 
High resolution X-ray diffraction (XRD) curves were then 
collected around the [004] symmetric Bragg reflection of 
silicon and are displayed in Fig. 2. A monochromatized 
CuK a i excitation source led to an angular resolution bet- 
ter than 5xl0 -3 degree. 

Data shown in Figures 1 and 2 illustrate the significant 
improvement of the structural quality and sample homo- 
geneity with respect to previous publications jTJ [9] . Up 
to 200 shots (cb < 4.5 at.%), we observe relatively flat 
boron concentrations (Fig. 1) with a sharp interface, at 
80 to 90 nm, on the order of 9 nm/decade. The peaks in 
XRD (Fig. 2) remain very well defined with a full width 
at half maximum on the order of 0.1°. As the number of 
laser shots is increased, the shift to higher diffraction an- 



FIG. 2: High-resolution XRD measurements of Si:B films at 
[004] Bragg reflection. The peak at 34.565° corresponds to 
the Si substrate and the second peak is characteristic of the 
doped layer for the indicated numbers of laser shots in GILD 
process. Inset: Relative variation of the lattice parameter 
perpendicular to the surface in % as a function of the boron 
concentration cb- Broken line, linear dependence (Vegard's 
law) at low doping. 



gles illustrates the contraction of the out-of-plane lattice 
parameter a,± as a result of the presence of boron atoms 
with a smaller covalent radius than silicon. This is a 
good indication of a sufficiently rapid re-crystallization 
limiting the diffusion processes and leading to an homo- 
geneous incorporation of boron on substitutional sites far 
above the solubility limit of 1.2 at%. In this range, as il- 
lustrated by the inset of Fig. 2, the out-of-plane negative 
strain deduced from the diffraction patterns varied lin- 
early with boron incorporation as measured by SIMS, the 
slope being in fair agreement with the strain rate coeffi- 
cient of Vegard's law in silicon [T3] . Above 200 shots, the 
interface began to exhibit a shoulder or even a large peak 
(500 shots) indicating a pileup of boron at the interface 
during the process. The diffraction peaks did not shift 
upwards anymore and remained close to 36° . Their inten- 
sity weakened and their width increased substantially, a 
second weaker broad peak even developing around 35.5° 
for 500 laser shots. This corresponds to a saturation of 
the out-of-plane strain around -3.5%. Since according to 
SIMS measurements the intake of boron went on increas- 
ing between 200 and 500 shots, one is led to conclude that 
in this range an increasing fraction of these boron atoms 
were not incorporated on substitutional sites. In our ex- 
perience, the saturation value of 36.1° observed here for 
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FIG. 3: Dependence of the superconducting transition tem- 
perature T c on the boron concentration cb deduced from 
SIMS measurements. In the inset : Resistive (right scale) 
and magnetic (left scale) transition in zero magnetic field in 
a Si:B film (400 laser shots and cb ~ 8 at.%) showing that 
both transport and diamagnetic transitions are very sharp 
and occur at the same temperature. 



FIG. 4: Temperature dependence of the upper critical field 
H C 2 deduced from transport measurements (see inset) in a 
heavily doped Si:B samples (c_b ~ 8 at.%, open squares; clas- 
sical theory in the very dirty limit, full line) together with 
the values previously reported in [I] (dots and dashed line). 
In the inset: temperature dependence of the resistance nor- 
malized to its normal state value in a serie of magnetic fields 
given in Gauss 



the [004] diffraction peak of Si:B is a quite general upper 
limit, independent of our GILD operation conditions. 

The superconducting temperatures and upper critical 
field have been deduced from ac-susceptibility (xac) and 
magneto-transport (R) measurements. In the Xac mea- 
surements, the films have been placed on top of miniature 
coils and T c has been obtained by detecting the change in 
the self induction of the coils induced by the supercon- 
ducting transitions. Both AC and DC magnetic fields 
were applied perpendicularly to the doped layer. Ap- 
proximate corrections were made to account for the de- 
magnetization factor but they don't affect the determi- 
nation of H C 2 where the magnetization is zero. As shown 
in the inset of Fig. 3, the onset of the diamagnetic re- 
sponse well coincides with the resistive transition, which 
is remarkably sharp in this sample, AT C ~ lOmK, while 
being only 5 times larger in the other samples of the se- 
ries. This allows an accurate determination of T c and the 
corresponding values have been reported in Fig. 3 as a 
function of the B content deduced from SIMS measure- 
ments. 

No superconductivity was found down to 40mK for the 



two less doped samples with c B ~ 1.4 at.% (50 laser 
shots) and Cb ~ 2.2 at.% (75 laser shots). For higher 
boron concentration, the critical temperature increases 
sharply with doping above some critical concentration c c 
~ 2.5 at.%o reaching 0.6K for eg ~ 8 at.%o. For compari- 
son, the sample measured in [T] displayed a roughly twice 
smaller T c with a much broader transition. The apparent 
saturation of T c is probably non-intrinsic since the incor- 
porated boron in this range of concentration might be 
non-elcctrically active as suggested by XRD and by the 
observation of a maximum angle of diffraction with broad 
and less intense peaks (see above). Calculations within 
Density Functional Theory and Virtual Crystal Approxi- 
mation [T3] predict T c ~ 0.3K (resp. 3K) for cb = 5 at.%o 
(resp. 10 at.%) using MacMillan formula [15]. Even if 
the order of magnitude is correct, the observed depen- 
dence of T c on cb is strikingly different from these pre- 
dictions. T c (cb) is actually very similar to that observed 
in doped diamond, with a critical threshold for doping 
concentration c c above which T c increases very sharply 
with a quasi power law dependence, T c cx (^— — 1) . 
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Whereas in doped diamond c c coincides with the critical 
concentration of the metal-insulator transition, in doped 
silicon c c ~ 2.5 at.%, i.e. about 300 times larger than 
the concentration necessary to induce metallicity. The 
meaning of c c in this latter case remains therefore rather 
puzzling and needs to be clarified by further studies. 

Electrical measurements were performed in magnetic 
field on the sample with cb ~ 8 at.% (400 laser shots) 
which has the highest T c and the sharpest transition in 
zero field. The results are displayed in the inset of Fig. 
4 as temperature was swept at different fixed magnetic 
fields. The transitions was shifted towards lower temper- 
ature as the magnetic field was increased, while remain- 
ing sharp and well defined. No hysteresis nor supercool- 
ing was observed, indicating a second-order transition, 
which is consistent with doped silicon being a type II 
superconductor. The upper critical magnetic field H C 2, 
defined by the usual criterion as where the resistance is 
90% of its normal state value, is plotted in Fig. 4. Note 
that the discussion below and the shape of the critical line 
don't depend on the exact criterion used; a R/R n =50% 
definition would only shift slightly the H c2 (T) line. 

The temperature dependence of H C 2 is very well re- 
produced within the standard microscopic theory. The 
full line in Fig. 4 is a fit to the solution of the linearized 
Gor'kov equations neglecting spin effects [TB] in the very 
dirty limit, | — ► 0, £ being the electronic mean-free path 
and £ the superconducting coherence length. Both the 
strong curvature and a 4 times larger H C 2(T=QK) ob- 
served previously [1] were most probably a consequence 
of the film inhomogcncity. Since the Fermi velocity vp is 
a parameter which cannot be determined experimentally 
with any accuracy, the electronic mean-free path i and 
H c 2(0) were estimated in a simple two-band free elec- 
tron model [17 with light (subscript Ih) and heavy (sub- 
script hh) holes having an effective mass mf h = 0.16 m e 
(resp. mttj = 0.49 m e ), m e being the bare electron mass, 
and with a total carrier concentration of 4xl0 21 cm~ 3 . 
The residual resistivity in this sample was equal to p ~ 
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The light hole band is expected to play a minor role 
for the determination of the upper critical field due to 
its larger Fermi velocity and its lower partial density 
of states. In the dirty limit, H C 2 is therefore related 
to the electronic mean-free path £ and to the coherence 
length in the heavy hole band by: H c2 w 2 ^f| hh 
[151 IH)] ~ 900 Gauss with £ hh being renormalized by 
the electron-phonon coupling constant, X ep ~ 0.3 [T] and 

^hh = jrA(i+A — ) ~ 1200nm with A being the su- 
perconducting gap. Considering the uncertainties in the 
thickness of the film and in the geometry of the con- 
tacts in a van der Pauw configuration together with the 
crudeness of a free electron approximation, this estimate 
H C 2(0) ~ 900 Gauss is in good agreement with the ex- 
perimental data H c2 (0) ~ 1050 Gauss (see Fig. 4). Note 

that the London penetration depth = y ^j^^fc ~ 
~ 60nm < £ implying that doped Si should be intrinsi- 
cally a type I superconductor such as SiC:B [5D] but is 
turned into a type II system by strong impurity effects 
(k ~ 0.7Xl/£ > 1 in the dirty limit but K = Ai/^<lin 
the clean limit). 

To conclude, the GILD technique is proved to be a 
powerful technique to dope silicon in the alloying range 
where superconductivity occurs. The superconducting 
transitions are sharp and well defined both in resistiv- 
ity and magnetic susceptibility allowing the study of the 
variation of T c on the boron concentration. This varia- 
tion is in contradiction with a classical exponential de- 
pendence on superconducting parameters. Instead, T c 
increases with a quasi-power law form above a critical 
threshold c c in striking similarity to what is observed in 
doped diamond. Although in the diamond case c c was 
equal to the MIT critical concentration, in silicon it lies 
far above. The comparison of these two cases offers a 
new playground for studying superconductivity in cova- 
lent systems near a MIT. 
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